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We study an area preserving parabolic map which emerges 
from the Poincare map of a billiard particle inside an elon- 
gated triangle. We provide numerical evidence that the mo- 
tion is ergodic and mixing. Moreover, when considered on 
the cylinder, the motion appear to follow a gaussian diffusive 
process. 
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The investigation of the quantum manifestations of 
classical dynamical chaos has greatly improved our un- 
derstanding of the properties of quantum motion. Even 
though, besides some very special cases, the non linear 
terms prevent exact solution of the Schrodinger equa- 
tion, still important useful information can be obtained 
concerning statistical properties of eigenvalues and eigen- 
functions. An important discovery has been the phe- 
nomenon of quantum dynamical localization |2| which 
consists in the quantum suppression of deterministic clas- 
sical diffusive behaviour. This suppression takes place 
after a relaxation time scale tji which is denned as the 
density p of the operative eigenstates Q , namely of those 
states which enter the initial conditions and therefore de- 
termine the dynamics. For times t < t^, the quantum 
motion mimics the classical diffusive behaviour and re- 
laxation to statistical equilibrium takes place. The re- 
markable fact is that quantum "chaotic" motion is dy- 
namically stable as it was illustrated in || . This means 
that, unlike the exponentially unstable classical chaotic 
motion, in the quantum case errors in the initial con- 
ditions propagate only linearly in time. More precisely, 
besides the relaxation time scale tn, a second very im- 
portant time scale exists, the so called random time scale 
t r ~ In h, below which also the quantum motion is expo- 
nential unstable. However, as remarked in |Q t r <C tR 
and therefore the quantum diffusion and relaxation pro- 
cess takes place in the absence of exponentially instabil- 
ity, ft should be noticed that, even though the time scale 
t r is very short it diverges to infinity as H goes to zero and 
this ensures the transition to classical motion as required 
by the correspondence principle. 

Therefore, typical quantum systems exhibit a new type 
of relaxation for which we do not have yet a physical 
description. In terms of the classical ergodic hierarchy, 
quantum systems can be at most mixing. While expo- 
nential instability is sufficient for a meaningful statistical 
description, it is not known whether or not it is also nec- 



essary. Several questions remain unanswered, e.g. there 
is no general relation between the rate of exponential 
instability and the decay of correlations. Moreover, as 
shown in B, quantum systems provide examples which 
show that linear dynamical instability is not incompati- 
ble with exponential decay of Poincare recurrences. 

In a recent paper a physical example has been found 
p|, a billiard in a triangle, which has zero KS entropy (the 
instability is only linear in time) but which possesses the 
mixing property H. This characteristics makes systems 
of this type, good candidates for the discussion of the 
above mentioned problems. In the present paper, start- 
ing from the discrete bounce map for the billiard in a 
triangle, we derive an area preserving, parabolic, clas- 
sical map. In others words, the map is marginally sta- 
ble i.e. initially close orbits separate linearly with time. 
We will show that this map is mixing, with power law 
decay of correlations and exponential decay of Poincare 
recurrences, and has a peculiar property: absence of pe- 
riodic orbits. Moreover, when the map is considered on 
the cylinder, it exhibits normal diffusion with the corre- 
sponding Gaussian probability distribution. 

Let us consider the following discontinuous skew- 
translation on the torus, with symmetric coordinates 
(x,!,)eT 2 = Kl)x[-l,l), 

Vn+i = Vn + asgnx„ +/3 (mod 2), 

x n+1 = x n + Vn+i (mod 2), (I) 

where sgnx = ±1 is the sign of x. The map (0), which we 
will call "the triangle map" is a parabolic, piece-wise lin- 
ear, one-to-one (area preserving) map, det J = I , tr J = 2 
with J := d{y n +x,x n +i)/d{y n , x n ) = 1. It is known that 
(continuous) irrational skew-translations (the above ma" 
(|f|) with a — and irrational (3) are uniquely ergodic 
and never mixing S, in fact they are equivalent to in- 
terval exchange transformations. However, the triangle 
map may have more complicated dynamics and we show 
below that discontinuity may provide a mechanism to 
establish the mixing property. Non-invertible piece-wise 
linear 2d parabolic maps have been studied in Ref. ||. 

The triangle map is related to the Poincare map of 
the billiard inside the triangle with one angle being very 
small. Indeed, let us assume that the small angle of the 
billiard can be written as 7 = tt/M with some integer 
M ^$> 1. Then the billiard dynamics may be unfolded by 
means of reflections over the two long sides of the triangle 
into the dynamics inside a nearly-circular 1M. -sided poly- 
gon. Within relative accuracy of 1/M, the approximate 



Poincare map inside such polygon, relating two succes- 
sive collisions with the short sides of the triangle — the 
outer boundary of the polygon — reads 

Un+l = V n + 2( Un - [u n ] - /j(-l) [u " ] ), (2) 

U n +1 = U n - 2v n+1 

where ju n is the polar angle and 7«„ is the angle of 
incidence of the trajectory in the n— th collision. The 
symbol [x] is the nearest integer to x. The parameter 
/i controls the asymmetry between the other two angles 
■n, ( of the triangle, namely 77, C = 7r/2 — 7(| ± /j) and 
we assume that the triangle has all angles smaller than 
7r/2, i.e. |/i| < 2- As shown in [gj, the system is equiv- 
alent to the mechanical problem of three elastic point 
masses on a ring (here one particle being much lighter 
than the other two). It is interesting to note that in the 
scaled variables (u, v) the small parameter 7 scales out 
from the map and the limit 7 — > simply means that the 
range of variables u n £ [0, 271-/7), v n £ [— 7r /(27),7r/(27)) 
becomes the entire plane R 2 . The above map can be 
compactified onto a torus T 2 by considering one 'prim- 
itive cell' (u n (mod 2),v n (mod 1)). After transform- 
ing the coordinates as y n = 2(—l) n (u n + v n ) (mod 2), 
x„ = (— l)"u„ — I (mod 2), we obtain the discontinous 
skew translation of a torus p) with a = 4/j and (3 = 0. 

In the following we consider the general case of the tri- 
angle map with parameters a and (3 being two indepen- 
dent irrationals. The particular case (3 = will be briefly 
discussed at the end of the paper. We fix the parameters 
value a = (|(V5- 1) - e" 1 )/^, (3 = (|(\/5- l) + e - 1 )/2, 
although qualitatively identical results were obtained for 
other irrational parameter values. 




t/N 

FIG. 1. The deviation of the filling rate r(t) from the 
random model in log-normal scale for three different mesh 
sizes N = 10 4 (dotted), N = 10 5 (dashed), and N = 10 6 
(solid curve). 

As a first step we make a detailed and careful test 
of ergodicity of the triangle map. To this end, follow- 
ing |12], we discretize the phase space T 2 in a mesh of 
N = Ni x Ni cells and then measure the number of 
cells n(t) visited by a given orbit up to discrete time 
t. Computing the phase space averages (.) by averag- 
ing over many randomized initial conditions we compare 



the quantity r(t) = (n(t)/N) thus obtained with the cor- 
responding trm(£) for the random model in which each 
throw onto a mesh of N cells is completely random. As it 
is known, in the latter case, rRjyi(i) = 1— cxp(— t/JV). The 
result shown in fig.p] provides strong evidence of (fast) er- 
godicity (without any secondary time scales): namely the 
exploration rate r(t) of phase space for the triangle map 
approaches 1 as t — > oo and, for sufficiently fine mesh N, 
is arbitrarily close to the random model 7"rm(£)- 
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FIG. 2. The time auto-correlation 

C(t) = (cos(7rj/t) cos(7T2/o)} averaged over 2 ■ 10 B orbits of 
length 16384 with randomized initial conditions. The dashed 
line has slope —3/2. In (b) we show the corresponding spec- 
tral density. Note that peak at ui = n/2 indicates a strong 
component of period 4. 

Having established with reasonable confidence that the 
triangle map is ergodic, we now turn our attention to the 
mixing property. This amounts to show asymptotic de- 
cay of time-correlation functions of arbitrary L 2 observ- 
ables. The extensive numerical experiments we have per- 
formed, suggest that arbitrary time-correlation functions 
decay asymptotically with a power law (f{t)g(0)) oc t~ a 
with the value of the exponent a close to a = 3/2. In 
figj|a we show the decay of auto-correlations of a typi- 
cal observable / = cos(7ry). ]13j The property of mixing 
and the nature of decay of correlations are intimately re- 
lated to the spectral properties of the unitary evolution 
(Koopman) operator over L 2 space of observables over 
T 2 . The value a > 1 we have empirically found implies 
absolutely continuous spectrum. Performing the inverse 
Fourier transform of C{t): 



C(t) := (/(t)/(0)) = / d» f {w)e 



(3) 



one calculates the spectral density d^if(uj)/duj which 
should be non-singular and continuous but non-smooth 
and non-analytic function, according to the (power-law, 
a > 1) nature of decay of correlations. In fig.|p we show 
the spectral density d/i f(ui) / 'dco which is apparently con- 
tinuous but not continuously differentiable function. In 
fact, we suggest that the discontinuities of the derivative 
are dense in order to ensure the correlation decay with 
the power a which is between 1 and 2. 
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FIG. 3. The Poincare recurrence probabilities P(t) for 
three different subsets [0, 0.1] x [0, 0.1] (solid), [0, 0.2] x [0, 0.2] 
(dashed), and [0,0.4] x [0,0.4] (dotted curve). Thick curves 
give numerical data obtained by computing the return proba- 
bility to the subset A for a single orbit of length 3- 10 11 . Thin 
curves are theoretical estimates for fully random dynamics, 
Pi(t) — exp(— fit), where fj, is the relative Lebesgue measures 
of the above sets, namely fi — 1/400, 1/100, 1/25 respectively. 

A very efficient tool for investigating the statistical 
properties of dynamical systems is the study of Poincare 
recurrences, i.e. the probability P(t) for an orbit to stay 
outside a specific subset A C T 2 for a time longer than 
t. In fig.|| we plot the Poincare recurrence probability 
P(t) for the map (fil) and for several different subsets 
of the form A = [0,6) x [0,6). The result is quite un- 
expected. Indeed, for any sufficiently small set (small 
6) the return probability appear to decay exponentially 
P(t) oc exp(— Xt). Moreover, the exponent A is very close 
to the Lebesgue measure of the subset /j, = \A\, as in 
the case of the random model of completely stochas- 
tic dynamics for which -Prm(^) = exp(— /if). Therefore 
the triangle map which is characterized by linear sepa- 
ration of orbits, exhibits exponential decay of Poincare 
recurrences, typical of hyperbolic systems. Notice that 
in strongly chaotic systems with positive Lyapunov ex- 
ponents, the presence of a zero measure of marginally 
unstable orbits (e.g. bouncing balls in the Sinai billiard) 
leads to a power law decay of Poincare recurrences. The 
simultaneous presence in our model of a power law decay 
of correlations and exponential decay of Poincare recur- 
rences is a fact for which, so far, we have no explanations. 
Indeed, even if there are no general rigorous theorems, it 



has been conjecutured that correlations of dynamical ob- 
servables have the same decay as the integrated Poincare 
recurrences |Q, namely C(t) ~ Pi(t) := / t °° drP{r) for 
asymptotically long times t. This relation is obviously vi- 
olated in our model (pi) and this interesting point requires 
further investigations. 
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P(t)-p(0) 
FIG. 4. The normal diffusion of the map (hi) on a cylinder. 
In (a) we show averaged squared displacement of an average 
over 10 7 orbits of length 2 ■ 10 5 (solid line) compared with 
the straight line (dashed) with slope D — 1.654. In (b) we 
show the corresponding distribution of displacements at three 
different times (t = 20000, 60000, 200000, solid curves) which 
are in perfect agreement with the solutions of the diffusion 
equation (Gaussians with variance Dt, dotted curves). 

Our last step is the investigation of the diffusive prop- 
erties of the system. To this end we consider the triangle 
map on the cylinder (y £ (— oo, oo)). In order to take into 
account the constant drift of y n with 'velocity' (3 we find 
convenient to introduce a new integer variable p n G Z: 



Vn = 2/o + (3n + ap n , 



(4) 



which has, by definition, vanishing initial value po = 0, 
and then study the diffusive properties in the variable 
p n . Our numerical results shown in (figjj) provide clear 
numerical evidence for normal diffusive behaviour. In 
particular we obtained a very accurate linear increase of 
the second moment (notice the long integration time) 

((Pn+t ~ Pn) 2 ) = (P 2 t ) = Dt (5) 

with diffusion coefficient D rs 1.654. The almost perfect 



Gaussian distributions of pt — Po — Pt obtained at differ- 
ent times (fig||b) indicate that we are in the presence of a 
normal Gaussian process. Note that dynamics (D) can be 
rewritten in terms of a closed map on an integer lattice Z 2 
with an explicit 'time-dependence', namely rewrite also 
the variable x n in terms of the integer variable q n G Z 

TLiTL -\- 1 J 

x n = Xo + y n + f3 Vaq n (mod 2). (6) 

and the map (|lj) becomes equivalent to an integer system 

Pn+l = Pn + (-l^+^on-h/S^+a^-^ (?) 

<7»i+l = q?i + Pn+l, 

with fixed initial conditions p = q a = 0. Here, the 
original initial conditions x$,yo enter as parameters. 

We note the trivial but important fact that the trian- 
gle map possesses no periodic orbits when the parameter 
j3 is irrational and a and j3 are incomensurable. There- 
fore, the general argument of Ref. Q using parabolic pe- 
riodic orbits cannot be used to derive the 1/t 2 decay of 
Poincare recurrence probabilities. It has been verified 
numerically that the non-existence of periodic orbits is 
indeed responsible for exponential decay of Poincare re- 
currence probability: When we replaced irrational with 
a crude rational approximation we have obtained a very 
clean crossover from initial exponential decay exp(— /it) 
to an asymptotic power law P(t) ex 1/t 2 due to existence 
of (long) periodic orbits. Our map thus provides a quite 
pathological example from the point of view of semiclas- 
sical periodic orbit theory, hence we pose an interesting 
question: which classical structure underpins the spec- 
tral fluctuations of the quantization of the triangle map 
(M? (see also Ref. 15 for skew translations, a = 0). 

An interesting special case of the triangle map is (3 — 
which, as discussed above, describes the dynamics of an 
elongated triangle (g). Here two cases should be distin- 
guished: i) the parameter a (= 4/i) is rational a = 2k/l, 
with k, I € Z, then the dynamics is pseudointegrable 
and confined onto ^-'valued' invariant curves (y n — yo)l 
(mod 2) = 0. ii) the parameter a is irrational, then the 
dynamics has been found to be ergodic. However ergodic 
properties turn out to be weak (see also [ fL0| ) and the rate 
of ergodicity is very slow as opposed to the general case 
(3 y^ 0: It has been shown that the number of different 
values of coordinate y n taken by a single orbit up to the 
discrete time T, < n < T, grows extremely slowly, as 
ex InT. A similar property has been found for triangular 
billiards in which one angle is rationally related with it, 
e.g. right triangles IbyJa]. In addition, numerically com- 
puted correlation functions of (fil) with (3 — 0, such as 
(cos(nyt) cos(iryo)}, show perhaps a tendency to decay as 
power-laws but with a small exponent a around 0.1. It is 
fair to say that, in this case, it is difficult to judge defini- 
tively, based on numerical experiments, on the property 
of mixing even though it cannot be excluded. 

In this paper we have shown that a gaussian diffu- 
sive process and mixing behaviour can take place in a 



simple area preserving map without dynamical exponen- 
tial instability. One may argue that parabolic maps are 
non generic and therefore irrelevant for the description 
of physical systems. However, the results presented here 
show that a meaningfull statistical description is possi- 
ble without the strong property of exponential instabil- 
ity. Even if the model discussed here is non generic in 
the context of classical systems, it can describe the typ- 
ical mechanism of quantum relaxation. Therefore it can 
play an important role in understanding and describing 
the quantum chaotic motion in analogy to the one played 
by Arnold cat map for classical systems. 
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